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Abstract. Whilst direct observations of internal transport in accretion disks are not yet possible, measurement of 
the energy emitted from accreting astrophysical systems can provide useful information on the physical mecha- 
nisms at work. Here we examine the unbroken multi-year time variation of the total X-ray flux from three sources: 
Cygnus X-1, the microquasar GRS1915-I-105, and for comparison the nonaccreting Crab nebula. To complement 
previous analyses, we demonstrate that the application of advanced statistical methods to these observational 
time-series reveals important contrasts in the nature and scaling properties of the transport processes operating 
within these sources. We find the Crab signal resembles Gaussian noise; the Cygnus X-1 signal is a leptokurtic 
random walk whose self-similar properties persist on timescales up to three years; and the GRS1915-I-105 signal is 
similar to that from Cygnus X-1, but with self-similarity extending possibly to only a few days. This evidence of 
self-similarity provides a robust quantitative characterisation of anomalous transport occuring within the systems. 

Key words, accretion disk - methods: statistical - X-rays: individual (Crab, Cygnus X-1, GRS1915-I-105) 
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1. Introduction 

Deeper understanding of the transport mechanisms that 
operate within accretion disks is important for a broad 
range of X-ray emitting astrophysical objects. In the ab- 
sence of local measurements, the key questions are (1) to 
quantify the ways in which the unseen transport processes 
are anomalous as distinct from diffusive; and (2) to estab- 
lish how this may be determined remotely from observa- 
tions of global quantities such as the input and outflow 
of energy. This information can then be used to provide a 
constraint for turbulence/instability models of astrophys- 
ical accretion disks. 

Here we analyse the total X-ray flux over several years 
from two accreting astrophysical objects - Cygnus X-1 
and GRS1915-I-105, the stellar-mass analogues of disk-jet 
active galaxies powered by a massive black hole - and, for 
comparison, the nonaccreting Crab which is powered by 
a neutron star. The reasons for selecting these sources for 
statistical analysis are twofold. First, since 1996 February 
20 they have been observed continually for several years by 
the All-Sky Monitor (ASM) on board the RXTE satellite 



(Swank et al. 2001 ), providing large data sets of around 
thirty thousand points. This enables us to seek correla- 
tions over several orders of magnitude up to the longest 
accessible timescales. Second, the source luminosities are 
sufficiently high to neglect instrument thresholds, uncer- 
tainties, and other sources in the field of view. The raw 
data are held at the Goddard Space Flight Center (GSFC) 
and can be accessed via their websiteR. Each point repre- 
sents the total X-ray flux (measured by the number of 
counts during periods that last 90 seconds) in the range 
1.3-12.2 keV, and the breakdown into three energy bands 
(1.3-3, 3-5 and 5-12.2 keV) is also available. We have 
analysed the different channels and found them to be sim- 
ilar to the total flux. Sampling intervals between the 90- 
second X-ray counting periods are distributed with means 
of 93 minutes for the Crab, 77 minutes for Cygnus X-1, and 
96 minutes for GRS1915-I-105. To give an indication of the 
spread of sampling intervals, 90% of the intervals arc below 
186 minutes for the Crab, 187 minutes for Cygnus X-1, and 
194 minutes for GRS1915+105. The implications for our 
techniques are explained where appropriate below (and see 
Appendix) . Calibration is undertaken by the ASM/RXTE 
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Fig. 1. X-ray time-series for the Crab, 1996 Feb. - 2001 ^ig. 3. X-ray time-series for GRS1915+105, 1996 Feb. 
Mar; mean timestep 93 minutes. 2001 Mar; mean timestep 96 minutes. 
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Fig. 2. X-ray time-series for Cygnus X-1, 1996 Sep. - 1999 
Dec; mean timestep 77 minutes. 



team and the processed data are freely accessible on their 
website^ (Bradt et al. ^OOl] ). The RXTE counts are not di- 
rectly proportional to luminosity, so the exact relationship 
between these two quantities would have to be accounted 
for in any model. However, this will affect only the nature 
of the PDFs and not the temporal correlation in the data. 
The three raw X-ray time-series are plotted in Figs. Il|, 
andy. 

Previous studies have focused on spatial and spec- 
tral structure, quasi-periodic oscillations, and modelling 
of temporal variability: see Weisskopf et al. ( 2000|) and 
references therein for the Crab; Maccarone et al. ( 2000 ) 



and re ferences therein for C ygnus X-1; and Belloni et al. 
(2000) and Rao et al. (2000) for summaries of the spectral 



and temporal analyses of GRS1915+105. Dhawan et al. 
( 2000| ) have studied the jet of this source, and the impor- 
tance of microquasars in general is discussed by Mirabel 
and Rodriguez ( |1999D . Nayakshin et al. ( |2000D model the 
gross variability of GRS1915-M05 on all but the shortest 



timescales, and recent observations by Chaty et al. ( |2001| ) 
have been used to search for interactions with the sur- 
rounding interstellar medium. As a complement to these 
techniques, we here examine three key statistical measures 
(described further in Sect. 2) for each time-series. These 
are the probability density function (PDF), the growth of 
range (extent of statistical self-similarity), and the differ- 
enced form, which we then compare with the well-known 
signatures of Gaussian noise and random walks. In partic- 
ular, we test for signatures of long-range correlations that 
are the hallmarks of turbulent transport. Insofar as these 
techniques can be applied successfully to the observational 
signatures of anomalous transport in the present context, 
they may also be transferable to related questions in space 
and laboratory plasma physics. 

There is already much evidence that accretion disks are 
locales for turbulent transport and instabilities. In order to 
explain typical accretion rates observed in a range of disk 



types, the standard disk model (Shakura & Sunyaev 1973 



Pringle & Rees 1972) uses turbulent viscosity to produce 
an appropriate outward transport of angular momentum. 
Following extensive numerical simulations (reviewed by, 
for example, Gammie 1998), the source of this turbulence 
is now believed to be a magnetic shear instability (Balbus 
& Hawley 1991). In addition, 1// power spectra (inverse 
power-law frequency dependence with unspecified index) 
- suggestive of highly-correlated, possibly self-organised 
critical, behaviour - have been observed in diverse ac- 
cretion systems (reviewed by, for example, Dendy et al. 
1998 ). It is also well established that instabilities occur- 



http://xte.mit.edu/XTE/asmlc/ASM.html 



ing in accretion disks (Dubus et al. 2001 and references 
therein) can give rise to anomalous statistics. However, 
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these statistics will differ from those of fully-developed 
turbulence, and hence we seek to constrain such models 
using the differencing and rescaling technique described 
below. 

Self-similarity, non-Gaussianity and non-trivial tem- 
poral scaling together are strong indications of highly- 
correlated processes such as turbulence (Bohr et al. IQQSj ). 
We will show how trivial scaling of near-Gaussian fluctua- 
tions in the Crab X-ray signal - evidence of diffusive trans- 
port - contrasts with non-trivial scaling of non-Gaussian 
fluctuations in the X-ray signals from Cygnus X-1 and 
GRS1915-t-105. Whilst there may be other methods we 
could apply, we are confldent that the methods described 
here are sufficiently insensitive to the timing and count- 
ing errors in the data; many techniques a re unr eliable even 
when random errors are small (Sornette |2000 ). 



2. Techniques 

2.1. The Probability Density Function (PDF) 

The first step in our analysis of each data set is the con- 
struction of its PDF. The PDF P{y) of a variable Y is de- 
fined such that the probability that Y lies within a small 
interval dy centred OTiY = y, is equal to P{y)dy. P{y) is 
normalised so that 



P{y)dy = 1. 



Here, we use dy ~ {ymax — ymm)/100 and further nor- 
malise each PDF to its mean {y) and standard deviation 
a to enable comparison with theoretical distributions. We 
note that the PDF was used by Antar et al. ( |200l| ) to char- 
acterise turbulent fluctuations in tokamak edge plasmas, 
and by Bhavsar & Barrow ( ^001 ) to model the magnitude 
distribution of the brightest cluster galaxies. Bramwell et 
al. < \im% ^000| , |2001| ) discuss the PDFs of fluctuations 
in highly-correlated systems (see Sect. 0), and Burlaga 
(2001) presents a review of log-normal distributions in the 
turbulent solar wind. 



2.2. Growth of range 

Consider a self-similar function y (t) . The difference be- 
tween the maximum and minimum values of y during a 
time interval At defines its range for that interval, R (At). 
Since y is self-similar, the ensemble-averaged value of R 
will scale with At. We may write 



{R (At)) = cAt" 



(1) 



with c and H constants; H here defines the Hurst expo- 
nent (Hastings & Sugihara 1993] ). For data that are only 
approximately self-similar, we use the above relation to 
check their closeness to self-similarity and also to obtain 
an effective value for H as follows. By moving the window 
At one point at a time through the raw data, an array 
of R (At) values is created from which the mean (R) is 
found (reducing the effects of uneven sampling). This is 



repeated for a range of At within the length of the data 
set. A plot of log(-R(Ai)) against logAi will reveal any 
deviations from self-similarity whilst the slope will give 
the best estimate of H. We use linear regression to calcu- 
late the 95% confidence interval for H (see Figs. 0, and 

Trivially, a function that is exactly constant over time 
has H = 0. At the other extreme, H = 1 indicates a func- 
tion whose range increases linearly with time (for positive 
c in Eq. H^). Intermediate values of H are generated by 
fractal functions, random Gaussian noise {H w 0.2), and 
Gaussian random walks (whose next value in time is the 
sum of the previous value and a random Gaussian incre- 
ment; H w 0.5). The value of H does not uniquely estab- 
lish correlation, however; uncorrelated series may present 
significant probabilities of observing greater values as the 
timescale increases. Consequently, the growth of range can 
be rather insensitive as a measure of correlation. We can 
in principle define a measure of correlation /3 in terms 
of fractal exponents such as H {(3 = 2H + 1). Following 
Malamud & Turcotte (1999), /3 = for uncorrelated noise 
and (3 — 2 for a Gaussian random walk. However, the 
use of only one method to estimate an unknown expo- 
nent (and hence /?) is to be avoided (Schmittbuhl et al. 
1995). For example, comparing values of H obtained from 
the growth of range of synthesized series with known /3, 
Malamud & Turcotte ( |1999D find that H « 0.2 for /3 < 0.5, 
whereas the slopes of the power spectra are equal to — /3 
for < /3 < 2. We therefore obtain, for comparison, es- 
timates of the strengths of correlation from the slopes of 
power spectra with Thompson multi-tapering (Percival & 
Walden 1993) to reduce the variance. 

The Hurst exponent has been u sed to quantify solar 

1997|) , correlation in 



magnetic complexity (Adams et al. 
the dynamics of the upper photosphere (Hanslmeier et 
al. 2000| ), and persistence in solar activity (Lepreti et al. 
POOO and references therein). In conjunction with rescal- 
ing techniques, the Hurst exponent has also been used to 
quantify self-similarity and long-range correlations in tur- 
bulen t fluc tuations in magnetic fusion plasmas (Carreras 
et aL |l99^) . 



2.3. Differencing and rescaling 

Starting from the raw data y (t) , we first form a set of 
differenced series Z {t, r) for a range of values of the time- 
lag r: 



Z{t,T)^y{t)-y{t-T). 



(2) 



From these we calculate a set of PDFs for Z{t), one for 
each value of r, which we denote by P{Z,t) (not nor- 
malised) - see Fig. 0a. It is often found to be the case 
that the peaks P (0, r) scale approximately as t~", as in 
Fig. ^. Such scaling is, of course, sought using the values 
of P (0, t) because they are the most accurate points of 
the PDFs, having the most data points in their bins. The 
question then arises whether the differenced series Z {t, r) 
can all be derived from a single PDF that scales. That is. 
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whether a common functional form emerges if we use the 
exponent a to rescale the axes such that 



Zt" = Z, 



and P -^ Pt- 



P. 



(3) 



causing the separate PDFs to collapse onto the same 
(r = 1) curve ~ Fig. Er (compare Mantegna & Stanley 
1995| ). Clearly, a reasonable straight line in Fig. Hb (at 
least up to some maximum timescale Tmax) is necessary 
but not sufficient for rescaling to succeed, since only the 
peaks of P {Z, r) are used to determine r. Physically, suc- 
cessful rescaling of the differenced X-ray data would im- 
ply, on all timescales up to Tmax , that (i) the sizes of X-ray 
fluctuations (that is, the differences between the observed 
values) are governed by a single type of process, and (ii) 
the total X-ray output is correlated, not random, in time. 
In the trivial case where P (Z, t) is independent of r, we 
would infer the absence of temporal correlation in the X- 
ray output. 

The value of a is given by a = —1/m where m is the 
slope of logP{0) v. logr as shown in Fig. 4(b). a char- 
acterises the common functional form of the distributions 
P(Z) viz: 

1. < a < 2 Levy (power-law with a -^ oo) 

2. a — 2 Gaussian (finite a) 

3. a > 2 power-law with finite a 



where P{Z) - |Z|-(i+") for Z -^ ±oo (Sornette ^OOOl ). 

Thus the differencing and rescaling procedure not only 
reveals the timescales over which physical processes oc- 
cur, but can also confirm any correlation suggested by 
the growth of range and inverse power-law form of the 
power spectra. Moreover, it quantifies the asymptotic be- 
haviour of the distribution of fluctuations, which is essen- 
tial for constraining turbulence/instability models (Bohr 
et al. |1998[). 



Differencing was used by Mantegna & Stanley (1995) 
to investigate fluctuations in the value of a financial index, 
but neither this technique nor the growth of range has yet 
(to our knowledge) been applied to astrophysical X-ray 
sources. 



o 



-1 



a. 
Cl -2 



-3 



i; 




-5 5 

(y-<y))/0 

Fig. 5. PDF for Crab X-ray time-series, normalised with 
respect to (y) and cr; dashed line is Gaussian. 
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Fig. 6. Logarithmic plot of growth of range R{At) (Eq. |l]) 
with At for Crab X-ray time- series; slope gives Hurst 
exponent 0.19 < H < 0.21. 



3. Results 

3.1. The Crab 

Fig. g shows the PDF of the total X-ray count-rate y (t) 
that form the time-series for the Crab (Fig. |^). The curve 
is close to Gaussian but with longer tails. Growth of the 
range is as low as that of Gaussian noise {H « 0.2, Fig. O); 
the slope of the power spectrum is better at detecting weak 
correlation (see Sect. |2.2[ ) and gives 0.26 < /3 < 0.28 at 
the 95% confidence level. Interestingly, the PDFs of the 
differenced data (Fig. m require no rescaling. That is, not 
only is the type of distribution independent of r, but the 
differences themselves show no spread over time. We infer 
that the raw time-series is uncorrelated; if the correlation 
suggested by the power spectrum exists, it is too weak to 
be detected by this method. 



3.2. CygnusX-1 

Fig. II shows the PDF of the raw count-rates from Cygnus 
X-1 (Fig. 0) over a continuous, relatively quiescent three- 
year period. In contrast to Fig. H, this PDF is clearly 
non-Gaussian. It appears possible to characterise some as- 
pects of its non-Gaussian statistical properties in terms of, 
for example, the Gumbel and Frechet distributions whose 
properties we outline in Sect. Meanwhile, we note from 
Fig. @ that the distribution of small amplitude events ap- 
pears to fit the left-hand tail of Gumbel distributions with 
a — [1,2] while the distribution of large amplitude events 
appears to fit the right-hand tail of a Frechet distribution 
with a — 1.1. This suggests that the total fiux has contri- 
butions from different physical processes, perhaps arising 
in different parts of the accretion disk and its surround- 
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Fig. 4. (a) Unsealed PDFs of differenced series (Eq. g), drawn from a Gaussian random walk y{t), each curve for a 
different t. (b) Power-law scaling of P(0,t) with r; -1/slope yields exponent a. (c) Rescaled PDFs (Eq. ^) share the 
same (Gaussian) curve for all t. 




Z=y(t)-y(t-T) 



Fig. 7. Unsealed PDFs of differenced time-series (Eq. ||) 
for Crab, with parameter r stepping up in half-integer 
powers of timesteps, to maximum 10^; dashed line is 
Gaussian. 



ings, whose individual PDFs could differ. Fitting to each 
of the many component distributions is problematic; here 
we simply note that the PDF in Fig. g could not be re- 
produced by summing Gaussian PDFs. 

Cygnus X-1 has a substantially higher growth of range 
{H w 0.35, Fig. |) than the Crab, but still below that of a 
Gaussian random walk {H « 0.5). This is confirmed by the 
power spectrum whose slope gives 1.00 < /? < 1.02 at the 
95% confidence level (see Sect. UM- This higher growth is 
evident in the unsealed PDFs of the differenced Cygnus 
X-1 data in Fig. [lO|. B y deriving a scaling exponent a from 
the peaks of Fig. [lO|, and using it to rescale as in Fig. |ll|, 
it is clear that the increments scale remarkably well over 
the full range of the data (four decades in time). This 
establishes both the existence of correlation in the X-ray 



(y-(y»/a 

Fig. 8. PDF for Cygnus X-1 X-ray time-series, normalised 
with respect to {y) and a; black dashed line Gaussian, 
green solid line log-normal, red solid line Frechet a = 1.1, 
blue dashed line Gumbel a — 1, blue solid line Gumbel 
a = 2. 



output, and variations controlled by one type of process, 
on timescales up to three years. However, with its slightly 
higher peak and broader tails, the PDF is distinctly non- 
Gaussian. In fact, the raw X-ray time-series of Cygnus 
X-1 is correctly described as a weakly leptokurtic random 
walk; that is, the P DF of the increments is long-tailed 
(Bouchard & Potters 2000| ). The asymptotic form of the 
increments P{Z) is quantified by a; in this case, 5.4 < 
a < 9.9 with 95% confidence. 
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Fig. 9. Logarithmic plot of growth of range R{At) (Eq. |l]) 
with At for Cygnus X-1 X-ray time-series; slope gives 
Hurst exponent 0.34 < H < 0.36. 




Z=y(t)-y(t-x) 

Fig. 10. Unsealed PDFs of differenced time-series (Eq. H) 
for Cygnus X-1, with parameter r stepping up in half- 
integer powers of timesteps, to maximum 10''. Curves with 
lower P(0) and broader tails correspond to higher values 
of r. 



3.3. GRS1915+105 

Fig. O shows the PDF of the raw data from 
GRS1915-f 105 (Fig. |). Like Cygnus X-1, GRS1915+105 
is better described in terms of non-Gaussian statistics than 
Gaussian. For example, the extreme left- and right-hand 
tails fit Frechet and Gumbel curves respectively, the re- 
verse situation of Cygnus X-1. Note also that there are 
several peaks, suggesting a multi-component source in- 
volving at least two physical processes with different sta- 
tistical properties. With H Ki 0.3 (Fig. |3|), midway be- 
tween Gaussian noise (H w 0.2) and a Gaussian ran- 
dom walk {H « 0.5), the growth of range of GRS1915's 
count-rate lies between those of the Crab and Cygnus 




Fig. 11. Rescaled PDFs of differenced time-series for 
Cygnus X-1, with parameter r stepping up in half-integer 
powers of timesteps, to maximum 10**; dashed line is 
Gaussian. Here, following Eq. Kl Zg denotes the differ- 
ences rescaled according to r and a, where a is obtained 
from values of P(0, r) in Fig. nO (compare Figs. Ha and 
ib); 5.4< a< 9.9. 
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(y-<y))/a 

Fig. 12. PDF for GRS1915-I-105 X-ray time-series, nor- 
malised with respect to (y) and a; black dashed line 
Gaussian, green solid line log-normal, red solid line Frechet 
a = 1.1, blue dashed line Gumbel a ~ 1, blue solid line 
Gumbel a = 2. 



X-1. This is confirmed by the slope of the power spec- 
trum: 0.91 < /3 < 0.93 at the 95% confidence level (see 
Sect. 2^). The PDFs of the increments do not rescale 
over the full four decades, see Fig. ^. The separate PDFs 
rescale approximately onto the same curve for r up to only 
1.5 decades (see Fig. |l^) with stronger leptokurtosis than 
for Cygnus X-1. However, we cannot unambiguously iden- 
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Fig. 13. Logarithmic plot of growth of range R{At) 
(Eq. |) with At for GRS1915+105 X-ray time-series; slope 
gives Hurst exponent 0.25 < H < 0.29. 
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Fig. 14. Rescaled PDFs of differenced X-ray time-series 
for GRS1915-f 105, with parameter r stepping up in half- 
integer powers of timesteps, to maximum 10^; dashed line 
is Gaussian. Here, following Eq. 0, Zs denotes the differ- 
ences rescaled according to t and a, where a is obtained 
from values of P(0, r) (compare Figs. Ha and 0b). Curves 
with lower -Ps(O) and broader tails correspond to higher 
values of r. 

tify scaling in this regime given the variation in sampling 
interval; see Appendix. 



4. Non-Gaussian statistical properties 

The tails of the PDFs of Cygnus X-1 and GRS1915-I-105 
are clearly non-Gaussian. Here, we compare these PDFs 
with those arising from extremal statistics. The two 
limiting distributions of interest are "Gu mbel's asymp- 
tote" an d Frec het (Fisher & Tippett |1928| ; Gumbel |l958 ; 
Sornette 200C). In outline, the limiting distributions re- 
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Fig. 15. Rescaled PDFs of differenced GRS1915+105 X- 
ray time-series, with r stepping up in single integer num- 
bers of timesteps, to a maximum value 10^'^; dashed line 
is Gaussian. 5.5 < a < 6.3. 



suit from selecting the maximum value ymax from each of 
a large number of large samples whose individual members 
are drawn from a distribution P{y). When P{y) decreases 
more rapidly than any power-law (as y — > oo), "Gumbel's 
asymptote" has the form 



Pciy^ax) = K{e^~'"T with 



u 



b{y-s) 



(4) 



where in the limit of an infinite number of measurements 
a = 1; the constants K, b, and s are fixed by normalisation 



as in Sect. 2.1 



Frechet distributions Ppiymax) arise in the same 
manner when the underlying PDF P(jj) is power-law. 
Mathematically, Ppiymax) is defined by Eq. ^ but with 
u = a + /31n(l -t- y/G), where K, h, s, a, and G are again 
fixed by normalisation as in Sect. 2.1, and /? = (1 — a)~^. 
These curves exist for 1 < a < 3/2. 

Physically, the fact that the tails of the Cygnus X- 
1 and GRS1915-I-105 data can be fitted to Gumbel and 
Frechet distributions may suggest that the observed sig- 
nals have the character of maximal events. In this case, 
they would be the brightest among multiple events (whose 
PDF could be power-law) occuring within each observa- 
tional time window. Interestingly, extremal statistics in a 
global measure are found in turbulent fluids and other 



highly-correlated systems (Bramwell et al. 1998, 2000, 
pool ). In these cases, anomalous values of a are found 
for Gumbel's asymptote, and we plot a = 1 and 2 along 
with the Frechet curve a — 1.1 for comparison with the 
X-ray data. Also, since our global quantity is emitted flux 
rather than absorbed power, our curves show the op posite 
handedness to the results of Bramwell et al. ( 1998[ ). 

We also plot the log-normal PDFs having the same (y) 
and (T as the data, and these curves fit as closely as ex- 
tremal distributions (see Figs. § and |l^). The significance 
and origin of extremal and log-normal distributions is cur- 
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rently of considerable interest in statistical physics and 
turbulence studies (Bramwell et al. 2001; Burlaga 2001). 



Appendix: Confidence limits for r 



5. Conclusions 



Using three key statistical methods - the PDF, the growth 
of range, and differencing and rescaling - we have iden- 
tified and quantified some fundamental contrasts in the 
character of the total X-ray output of the Crab, Cygnus X- 
1, and GRS1915-f 105 over multi-year time intervals. The 
Crab shows near-Gaussian behaviour with its low Hurst 
exponent and lack of temporal scaling. This is to be ex- 
pected since the flux emanates from a very wide region 
over which correlation is unlikely, there being no obser- 
vational evidence of an accretion disk in this source. In 
contrast, Cygnus X-1 has a time-series with scaling over 
three years, and resembles a random walk with leptokur- 
tic increments whose PDF satisfies P{Z) ~ |Z|~'^^+"^ for 
Z -^ ±oo and a « 8. Its PDF closely follows multi- 
component Gumbel and Frechet curves, suggestive of the 
dominance of maximal events. Similarly, the output from 
GRS1915-I-105 lies close to such curves and its increments 
are also leptokurtic, although the roles of Frechet and 
Gumbel components in the tails are reversed. However, 
as suggested by a lower Hurst exponent, results from this 
source are consistent with only short-range scaling with 
a PS 6, and more data are required to establish this find- 
ing. Thus we have evidence that the two accreting objects 
display a degree of correlation in their X-ray time-series, 
which is absent from the nonaccreting Crab. This is a 
quantitative, observational, and model-independent mea- 
sure of anomalous (non-diffusive) transport in accretion 
disks. 
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Table 1. Median-centred 90% confidence limits for 
log 10(t) in the X-ray time series of the Crab, Cygnus X- 
1, and GRS1915-I-105. When considering fluctuations over 
small numbers of data points, a consequence of the uneven 
observational sampling times is that a broad range of ac- 
tual times are represented by each value of r (see Fig. [l5|) . 
To quantify these ranges, this table lists the limits between 
which 90% of the actual times lie for each value of r. 
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